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A  finite  element  technique  is  presented  for  the  plastic  analysis  of 
structures  subjected  to  out-of-plane  bending  alone,  or  in  combination 
with  in-plane  membrane  stresses.  The  method  makes  use  of  a  linear 
matrix  equation  of  finite  element  analysis,  formulated  to  include  the 
effect  of  initial  strains.  This  equation  is  applied  to  the  plasticity 
problem  by  Interpreting  plastic  strains  as  initial  strains,  the  material 
nonlinearity  being  introduced  throu^  subsidiary  stress-strain  re¬ 
lations  from  an  Incremental  plasticity  theory.  In  addition,  the  analysis 
is  combined  with  an  incremental  technique  developed  to  account  for 
the  effects  of  geometric  nonlinear  behavior.  Thus,  the  present  analysis 
is  capable  of  treating  the  combined  effects  of  material  and  g^metric 
nonlinearity.  Application  of  the  procedure  is  made  to  beam  and  arch 
structures  in  the  presence  of  both  types  of  nonlinearity,  and  to 
rectangular  plates  for  which  material  nonlinearity  alone  is  present. 
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SECTION  I 
INTRODUCTION 


The  description  of  plastic  behavior  presents  some  basic  difficulties  to  the  structural 
analyst.  These  difficulties  are associatedwithaproper description  of  the  material  phenomenon 
and  the  nonlinear  nature  of  the  resulting  governing  equation.  Thus,  the  mathematical  formu¬ 
lation  of  the  plasticity  problem  makes  an  analysis  of  aU  but  the  simplest  structures  a  very 
formidable,  if  not  an  Impossible,  task.  Consequently,  considerable  attention  has  been  given 
recently  to  the  extension  of  finite  element  techniques  to  Include  the  effects  of  plastic  behavior 
(References  1-5).  These  techniques  have  the  advantage  of  being  capable  of  treating  the  effects 
of  plasticity  in  complex  structures  by  .utilizing  various  algorithms  for  linearizing  the  basic 
nonlinear  nature  of  the  problem. 

Most  of  the  current  effort  concerned  with  the  application  of  finite  element  techniques  to 
the  plasticity  problem  has  been  limited  to  the  treatment  of  structures  in  states  of  membrane 
stress.  In  addition,  the  methods  generally  neglect  the  effects  of  geometric  nonlinearity.  These 
limitations  are  too  restrictive  for  many  Important  aerospace  structures.  Consequently,  it  is 
the  purpose  of  the  present  paper  to  extend  the  methods  already  developed  to  provide  for  a 
plastic  bending  analysis  that  accounts  for  membrane  stress  states  and  geometric  non¬ 
linearity. 

The  present  method  makes  use  of  a  governing  linear  matrix  equation  that  relates  the 
applied  loading  to  the  nodal  displacement  and  initial  strains.  For  the  purpose  of  a  plas¬ 
ticity  analysis,  the  plastic  strains  are  interpreted  as  initial  strains.  Use  of  the  initial  strain 
concept,  to  treat  the  effects  of  plasticity,  requires  the  development  of  appropriate  matrix 
relations  based  on  assumptions  for  the  distribution  of  both  displacements  and  initial  (plastic) 
strain  within  a  finite  element.  The  specification  of  a  distribution  for  plastic  strain  within  a 
finite  element  forms  the  basis  on  which  the  present  plasticity  analysis  depends. 

Inclusion  of  the  effects  of  geometric  nonlinearity  is  primarily  of  concern  in  problems 
involving  thin  beams,  plates,  and  shells  in  the  plastic,  as  well  as  elastic,  range.  A  finite 
element  method  that  utilizes  an  incremental  procedure  requiring  a  successive  modification 
of  the  element  stif&iess  properties  has  been  discussed  in  Reference  6.  This  method  requires 
the  introduction  of  an  additional  stiffness  matrix  to  account  for  the  effects  on  the  bending 
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stiffness  of  the  membrane  stresses  generated  as  a  consequence  of  geometric  nonlinearity. 
In  addition,  the  effect  that  changes  in  geometry  have  on  subsequent  deformations  is  taken  Into 
account.  This  Incremental  procedure  is  incorporated  Into  the  plastic  bending  analysis  to 
treat  the  combined  effects  of  material  and  geometric  nonlinearity. 


SECTION  II 

MATERIAL  NONLINEARITY-PLASTIC  BENDING  ANALYSIS 

An  important  advantage  of  finite  element  techniques  is  the  ability  to  specify  the  distri¬ 
bution  of  displacement  and  strain  states  within  each  finite  element.  This  permits  assumptions 
to  be  made  for  the  distribution  of  plastic  strain  and  the  development  of  regions  of  plasticity 
within  an  element.  These  assumptions  considerably  reduce  the  complexity  of  the  analysis  by 
defining  the  distribution  of  plastic  strain  in  any  element,  one  the  nodal  values  are  deter¬ 
mined.  This  feature  is  consistent  with  finite  element  analysis,  and  allows  us  to  be  concerned 
only  with  quantities  at  node  points  of  the  idealized  structures. 

For  the  case  of  membrane  stress  states,  the  plastic  strains  are  assumed  to  vary  in  a 
prescribed  manner  in  the  plane  of  the  element.  For  out-of -plane  bending,  an  assumption  must 
be  made  for  the  distribution  of  plastic  strain  throu^  the  thickness  as  well  as  in  the  middle 
surface  of  the  element.  Specifically,  the  present  analysis  assximes  the  plastic  strains  to  vary 
linearly  along  the  edges  of  a  finite  element  between  adjacent  nodes,  and  in  addition  assumes 
a  linear  variation  of  plastic  strains  from  the  upper  or  lower  surface  of  the  element  to  an 
elastic -plastic  boundary  (or  boundaries)  located  within  the  cross  section  of  the  element. 
These  assumptions  require  the  determination  of  the  position  of  an  elastic-plastic  boundary 
based  on  its  assumed  distribution  within  the  element  during  the  course  of  loading.  Thus,  the 
present  analysis  utilizes  the  concept  of  a  finite  element  in  which  there  is  a  progressive 
development  of  a  plastic  region  Instead  of  the  layered  approach  of  Reference  7  or  the  sandwich 
idealization  of  Reference  8. 

The  above  assumptions  are  made  in  the  development  of  the  governing  linear  matrix 
relation,  which  has  been  formulated  to  Include  the  effects  of  Initial  strains.  These  assumptions, 
as  applied  to  a  typical  beam  finite  element  for  which  pure  bending  behavior  has  been  assumed. 
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are  shown  In  Figure  1.  The  function  for  the  displacement  in  the  z-dlrecUon  Is  assumed  to  be 

of  cubic  order  in  the  coordinate  x,  and  is  written  in  terms  of  the  generalized  nodal  dls- 
placements  as 


In  choosing  a  displacement  function,  it  is  important  to  include  all  fundamental  strained  states 
and  all  rigid  body  terms.  Equation  1  satisfies  these  requirements  for  a  beam  element,  and 
allows  for  a  constant  shear  load  and  linearly  varying  moment  along  the  length  of  the  element. 
The  plastic  strain  distribution  Is  assumed  to  vary  linearly  in  the  x-dlrectlon  from  its  value 
at  the  upper  (or  lower)  surface  at  node  1,  represented  in  Figure  1  as  to  its  value  at  the 
upper  (or  lower)  surface  at  node  j,  represented  as  This  assumed  distribution  is  written 


where  z  represents  the  depth  of  the  elastic-plastic  boundary.  In  addition,  as  seen  from 
Equation  2,  it  is  assumed  that  at  a  node,  the  plastic  strain  varies  linearly  from  its  value  at 

the  upper  or  lower  surface  to  zero  at  an  elastic-plastic  boimdary  located  through  the  cross 
section. 


'Oi 


(2) 


The  depth  of  the  elastic-plastic  boundary,  which  propagates  from  the  upper  and  lower 
surface,  is  measured  from  the  neutral  axis  for  pure  bending,  as  shown  in  Figure  1.  In  general, 
the  depth  of  this  boundary  cannot  be  directly  related  to  the  load.  Hence,  the  value  of  I  must 
be  determined  from  the  total  strain  distribution,  which  is  assumed  to  vary  linearly  throu^ 
the  thickness  in  accordance  with  Klrchoff's  hypothesis.  The  functional  form  representing 
the  distribution  of  the  elastic-plastic  boundary  is  assumed  to  be  a  linear  function  of  the 
coordinate  x  and  may  be  written  as 

‘  13) 

where  z^  and  represent  the  depth  of  the  elastic -plastic  boundary  at  nodes  1  and  J,  re¬ 
spectively.  Thus,  with  the  preceding  assumptions,  the  elastic -plastic  boundary  consists  of 
a  surface  in  the  interior  of  the  element  that  extends  over  the  entire  area  of  the  element  and 
Intersects  the  edges  along  strai^t  lines  Joining  nodes,  as  Illustrated  in  Figure  1.  In  addition. 
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these  assumptions  eliminate  the  necessity  of  determining  an  elastic-plastic  boundary  on  the 
faces  of  the  element  between  nodes,  but  still  require  locating  such  a  boundary  throu^  the 
thickness. 

The  present  assumptions  are  further  extended  to  Include  the  effects  of  bending  in  com¬ 
bination  with  a  membrane  stress  state.  As  seen  In  Figure  2,  this  extension  necessitates  the 
determination  of  the  position  of  the  elastic-plastic  boundary  relative  to  both  the  iqjper  and/or 
lower  surface.  The  functional  representation  of  the  plastic  strain  distribution  and  the 
representation  of  the  elastic -plastic  boundary  are  taken  similar  to  Equations  2  and  3  written 
for  both  the  upper  and  lower  surface.  A  second  matrix,  in  addition  to  the  usual  stiffness 
matrix,  termed  the  initial  stress  stiffness  matrix  (to  be  discussed  in  the  next  section),  is 
also  implemented  to  account  for  the  effects  of  the  membrane  load  on  the  bending  stiffness. 
This  problem  also  requires  the  introduction  of  a  second  displacement  component,  u,  acting 
in  the  axial  direction, 


(4) 


It  should  be  noted  that  although  the  functional  form  of  the  plastic  strain  distribution,  as 
shown  in  Figure  2,  does  assume  the  existence  of  a  neutral  axis  located  within  the  cross  section 
of  the  beam  element,  the  present  analysis  is  capable  of  considering  plastic  sections  in  which 
the  neutral  axis  is  not  located  within  the  thickness  of  the  beam;  that  is,  the  upper  and  lower 
strains  are  of  the  same  sign.  This  situation  occurs  with  the  application,  or  generation,  of 
large  membrane  stresses  as  compared  to  the  existing  bending  stresses.  The  treatment  of 
this  situation  is  accomplished  by  modifying  the  functional  form  of  the  plastic  strain  distri¬ 
bution  given  in  Figure  2. 

The  present  method  has  also  been  extended  to  treat  the  more  complex  problem  of  the 
plastic  bending  of  a  plate.  A  typical  rectangular  plate  element  Is  shown  In  Figure  3.  The 
displacement  function  chosen  is  the  one  originally  used  by  Bogner,  Fox,  and  Schmit  (Ref¬ 
erence  9),  and  is  in  terms  of  products  of  first  order  Hermltian  pol3momials.  The  conq>onents 
of  initial  strain  are  assumed  to  vary  as  products  of  zero  order  Hermitlan  pol3momlals  In  the 
plane  of  the  element  and  linearly  through  the  cross  section  from  their  values  at  the  upper 
(or  lower)  surface  to  zero  at  the  elastic-plastic  boundary.  The  depth  of  the  boimdary  through 
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f 


Figure  1.  Typical  Elastic-Plastic  Beam  Element  (Pure  Bending^ 


Figure  2.  Typical  Elastic-Plastic  Beam  Element  (Bending  and  Axial  Loads) 
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Figure  3.  Typical  Elastic -Plastic  Rectangular  Plate  Element  (Pure  Bending) 
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the  thickness »  which  must  be  determined  at  each  of  the  four  nodes  of  the  rectangular  element* 
is  computed  from  the  total  strains  by  means  of  the  following  relation: 


and  the  superscript  T  denotes  total  strains.  The  function  defining  the  elastic -plastic  boundary 
in  the  plane  of  the  element  is  assumed  to  be  in  the  form  of  products  of  zero  order  Hermitian 
polynomials  as  shown  in  Figure  3.  The  foregoing  assumptions  associated  with  the  plastic 
strain  distribution  and  the  representation  of  the  elastic-plastic  boundary  ensure  compatibility 
of  these  quantities  along  element  boundaries. 

METHOD  OF  ANALYSIS 

Once  the  assumption  is  made  on  the  distribution  of  displacement*  the  total  strain  distri¬ 
bution  can  be  e^ressed  in  terms  of  nodal  displacements  by  making  use  of  the  appropriate 
strain  displacement  relations  in  conjunction  with  the  assumed  displacement  function.  This 
relation  can  be  written  in  matrix  form  as  follows: 

/  =  w  (5) 

0 

where  t  ^  is  the  vector  of  total  strains 

is  the  vector  of  generalized  nodal  displacements. 

The  assumed  distribution  of  plastic  strains  can  be  written  In  terms  of  their  nodal 
quantities  as 


€ 


where  is  the  vector  of  nodal  plastic  strains. 


(6) 
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The  matrices  W  and  are,  In  general,  functional  matrices  which  depend  explicitly  on 
the  assumptions  made  for  ^e  distribution  of  displacement  and  Initial  strain,  respectively, 
within  the  element.  Specifically,  the  assumptions  discussed  above  for  the  beam  and  plate 
elements  are  used  In  the  formation  of  W  and  Wp  in  the  present  analysis.  The  element  stiffness 
properties  in  the  presence  of  initial  strains  are  developed  on  the  basis  of  these  two  func¬ 
tional  matrices.  They  are  obtained  by  substituting  Equations  5  and  6  Into  the  e;qpresslon  for 
strain  energy  and  then  employing  Castlgllano*s  first  theorem.  To  this  end,  the  expression 
for  strain  energy  which  e3q)llcltly  contains  the  effect  of  initial  strains  may  be  written  as 

u  =  -1  E  J  dV-  E  <  dV  (7) 

V  V 

where  the  elements  of  the  matrix  E  are  the  coefficients  associated  with  the  appropriate 
linear  stress  strain  relations.  Substitution  of  Equation  5  and  the  assumed  plastic  strain 
distribution  of  Equation  6  into  Equation  7  leads  to 

u  =  Y  (Iq  k  do  -  <1q  K*  €o  ‘8) 

where 

k=///w’EWdV  (9) 

V 

l.*  =  ///WE  Wp  .V  1,0, 

''p 

The  matrix  k  represents  the  element  elastic  stiffness  matrix  and  depends  only  on  the 
assumption  made  for  element  displacements.  The  second  Integral,  Equation  10,  represents 
the  initial  strain  stiffness  matrix  which  Is  dependent  on  the  assumptions  for  both  total  and 
plastic  strains.  The  quantity  Vp  in  Equation  10  Is  the  volimie  of  the  plastic  region  In  each 
element  as  determined  by  the  representation  of  the  elastic -plastic  boundary.  Consequently, 
the  elements  of  the  Initial  strain  stiffness  matrix  k  *  will  be  a  function  of,  among  other 
quantities,  the  depth  of  the  elastic-plastic  boundary  at  each  node, and  must  therefore  be 
continuously  computed  during  the  course  of  loading.  The  Initial  strain  matrix  for  the  beam 
element  subjected  to  pure  bending  and  that  for  combined  bending  and  axial  load  are  given  In  the 
Appendix,  Due  to  space  limitations  the  initial  strain  matrix  associated  with  the  rectangular 
plate  bending  element  is  omitted  here.  It  may  be  found  In  Reference  14. 
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Application  of  Castigliano’s  theorem  to  Equation  7  yields  the  governing  linear  matrix 
equation  for  an  individual  finite  element 


dU 

““o 


(ID 


where  Pq  is  the  vector  of  generalized  nodal  forces. 


A  similar  equation  is  also  developed  In  incremental  form  in  anticipation  of  combining 
plasticity  with  the  incremental  geometric  nonlinear  analysis.  This  equation  is  written  as 
follows: 


AP-  =  k  Ad_  -  k  A  c_ 

0  0  0 


(12) 


It  should  be  emphasized  that  the  initial  strain  matrix  associated  with  the  increment  of 
plastic  strain  A  f  q  in  Equation  12  is  written  as  k  to  distinguish  it  from  the  inital  strain 
matrix  k  *  associated  with  the  total  plastic  strain  .  These  matrices  may  differ  sub¬ 
stantially,  since  the  functional  form  assumed  for  the  distribution  of  total  plastic  strain  will 
not,  in  general,  coincide  with  the  assumed  distribution  associated  with  the  Increment  of 
plastic  strain.  In  the  present  analysis,  was  assumed  to  vary  linearly  throu^  the  thickness 
from  the  upper  and/or  lower  surface  to  the  elastic-plastic  boundary.  This  assumption  implies 
a  bilinear  distribution  of  Acq  .  The  specification  of  a  distribution  of  this  form  requires  the 
determination  of  a  value  of  the  plastic  strain  increment  at  some  intermediate  point  in  the 
cross  section,  in  addition  to  a  value  at  the  upper  and/or  lower  surface. 


To  avoid  the  added  complexities  associated  with  the  use  of  Equation  12,  we  may  alter¬ 
natively  use  an  Incremental  form  of  Equation  11  as  follows: 


Af^  =  k  Ad*^  -  Aq 
0  0 

where 


03) 


Aq  =  A 


l-l 


o 


and  the  superscripts  i  and  i-1  refer  to  the  current  and  preceding  loading  step,  respectively. 
In  Equation  13,  the  product  of  the  Initial  strain  matrix  and  the  total  plastic  strain  is  con¬ 
sidered  as  a  vector  of  fictitious  loads.  The  Increment  of  this  fictitious  load,  represented  as 
Aq  •  is  determined  at  any  step  by  subtracting  its  current  value  from  that  computed  in  the 
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preceding  step.  In  this  manner,  only  total  values  of  plastic  strain  are  utilized  in  the  governing 
linear  matrix  equation.  The  desired  form  of  the  equation  is  obtained  by  grouping  the  in¬ 
crements  of  generalized  nodal  forces  and  fictitious  forces,  resulting  in  the  following  equation: 


k  Ad^'  =  AP'  +  Aq'"  ' 


(14) 


Here,  it  is  seen  that  the  values  of  the  Increment  of  fictitious  load  introduced  into  Equation  14 
are  values  taken  to  be  equal  to  those  computed  in  the  preceding  load  increment.  The  use  of 
this  type  of  predictor  procedure  is  necessary  because  the  depth  of  the  elastic-plastic  boxmdary 


(and  the  cxirrent  value  of  plastic  strain)  at  those  nodes  of  the  structure  in  the  plastic  range 


can  be  determined,  in  general,  only  from  the  stress  (or  strain)  distribution  computed  at  the 
end  of  the  load  step.  The  position  of  the  elastic -plastic  boundary  is  determined  at  the  end  of 
each  load  increment,  and  is  assumed  to  remain  fixed  dxudng  the  next  Increment. 

Equation  14  is  written  for  each  element  in  the  structural  idealization  and  then  appropri¬ 
ately  assembled  to  form  the  overall  linear  matrix  equation  for  the  entire  structure.  This 
equation,  not  shown  here,  is  identical  in  form  to  Equation  14.  The  incremental  solution  tech¬ 
nique  using  this  equation  reduces  to  a  sequence  of  linear  problems  in  which  the  applied  loading 
is  constantly  modified  by  the  fictitious  force  vector, 

PLASTICITY  RELATIONS 

The  foregoing  matrix  equations  must  be  used  in  conJ\mction  with  an  appropriate  plasticity 
theory.  Plasticity  is  introduced  into  Equation  14  through  subsidiary  stress-strain  relations 
provided  by  this  theory.  In  the  present  paper,  consideration  is  given  to  both  elastic,  strain- 
hardening  or  Ideally  plastic  material  behavior.  For  both  types  of  behavior,  the  total  strain 
increment  at  a  node  canbe  written  as  the  sum  of  an  elastic  and  plastic  con^)onent,  represented 
as  A«q  and  ,  respectively,  in  the  following  equation: 


A  el.  =  A«*  +  Af* 
0  0  0 


(15) 


In  addition,  the  increment  of  elastic  strain  is  related  to  the  stress  increment,  A  <r  ,  by  means 
of  Hookers  Law.  Thus,  Equation  15  can  be  written  as 


(16) 
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For  an  elastic,  strain-hardening  material,  we  make  use  of  a  linear  Incremental  relation 
between  plastic  strain  and  stress  increment 

Ac  =  C  Act  (I7) 

This  relationship  is  represented  in  a  general  form  by  the  matrix  c  in  Equation  17.  The 
formulation  of  this  matrix  is  directly  related  to  the  plasticity  theory  chosen  for  use,  l.e., 
these  elements  may  be  determined  by  using  an  isotropic  hardening  theory  or  the  kinematic 
hardening  theory  of  plasticity.  The  elements  of  this  array  for  plane  stress,  obtained  by 
using  Drucker’s  postulate  with  the  Prager -Ziegler  kinematic  hardening  theory,  are  explicitly 
given  in  Reference  1. 


Substituting  Equation  17  into  Equation  16  leads  to  an  Incremental  stress-strain  relation 
given  in  the  Equation  18 


Act  s  R  *  A«^ 


(18) 


where 

R  s  e”'  +  C 


It  should  be  noted  that  there  is  no  unique  stress  increment  for  a  given  plastic  strain  incre¬ 
ment  vector.  Therefore,  the  matrix  C  ,  given  in  Reference  1,  is  singular.  However,  the 
matrix  R  ,  defined  in  Equation  18,  will  possess  an  inverse,  thereby  providing  the  necessary 
coefficients  relating  the  stress  increment  to  the  increment  of  total  strain. 


The  increment  in  total  strain  at  anode,  A  e^,  is  obtained  from  the  increment  of  displace¬ 
ment  by  using  Equation  5  in  incremental  form  as  follows: 


*0 


09) 


where  Wq  Is  defined  at  a  node.  It  should  be  noted  from  Equations  19  and  5  that  the  functional 

form  chosen  to  represent  the  increment  of  total  strain,  A  •  ,  is  Identical  to  that  used  in  the 

T  ^ 

representation  of  its  full  value,  «q.  Thus,  having  obtained  the  increment  of  displacement 
from  the  solution  of  the  total  linear  matrix  equation  in  the  form  of  Equation  14,  Equations  15 
through  19  represent  the  necessary  relations  that  must  be  used  to  obtain  the  complete 
solution  for  increments  of  stress  and  strain,  assuming  elastic  strain-hardening  material 
behavior.  After  summing  all  incremental  quantities,  a  new  value  of  the  Increment  of  ficti¬ 
tious  load,  Aq  ,  is  determined  for  each  element  in  the  plastic  range  and  the  procedure  is 
repeated  until  the  desired  maximum  value  of  the  load  is  reached. 
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Consideration  of  elastic,  ideally  plastic  material  behavior  is  necessary  for  predicting 
the  collapse  load  of  a  given  structure  and  loading  situation.  The  two  conditions  to  be  satisfied 
for  multiaxial  elastic,  ideally  plastic  material  behavior  are: 

1)  the  stress  increment  vector  must  remain  tangent  to  the  loading  surface,  and 

2)  the  plastic  strain  Increment  vector  must  remain  normal  to  the  loading  surface, 
where  the  loading  function  is  the  representation,  in  stress  space,  of  the  initial  yield  function, 
or  the  subsequent  yield  function  after  some  plastic  deformation  has  occurred. 

The  above  conditions  are  expressed  analytically  for  the  case  of  plane  stress  and  using  the 
von  Mises  yield  condition  by  the  following  two  equations: 


(20) 


(21) 


where  d  X  is  a  positive  scalar  quantity. 

If  we  consider  the  differential  of  the  stress  components  as  Incremental  quantities,  the 
impliclty  differential  of  Equation  20  provides  a  linear  relation  among  the  components  of 
stress  increment,  represented  in  matrix  form  as  follows: 


Aor  =  E  Aor 


(22) 


Expressing  A  in  terms  of  A  c  ^  and  A  yields  the  elements  of  the  matrix  E  as 


0  -  M.  -M, 


2 


E  =  0 


0 


0  0 


where 


Mg  =  “  T  * 
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If  we  replace  dCy  by  Ac.j,  Equation  17  provides  a  linear  relation  between  components  of 
plastic  strain  Increment,  written  as 


A<  =  E  A« 

Pa  ^ 

Then  e^ressing  Ac  and  Ay  in  terms  of  Ac  gives  the  elements  of  the  matrix  E  , 

y  ^  * 


(23) 


E  = 


I 

M, 


0 

O 

0 


0 

0 

0 


It  is  apparent  from  Equations  20  and  21  that  only  two  of  the  three  components  of  stress  in¬ 
crement  and  only  one  of  the  three  components  of  plastic  strain  increment  are  required  to 
obtain  the  remaining  components.  Thus,  only  three  of  the  six  quantities  are  independent 
variables.  The  increments  of  stress  and  plastic  strain  can  now  be  written  in  terms  of  a 

vector,  A  <(>  ,  representing  these  Independent  quantities,  aiisitrarily  chosen  as  Ac  ,  Acr  , 

X  y 


and  At 


xy' 


A«  =  <  Ao; 

At, 


(24) 


xy 


Equations  22  and  23  may  now  be  rewritten  to  relate  the  increments  of  stress  and  plastic 
strain  In  terms  of  Aw 

Ac  =  E  A« 


and 


A<  =  Aw 


(26) 


A  relation  between  the  vector  A  cj  and  the  Increment  of  total  strain  is  obtained  by  sub¬ 
stituting  Equations  25  and  26  into  Equation  16. 


.  T  ♦  . 

At  =  E  Aw 


(27) 


where 


♦  “I  -  ^ 

E  5  E  E  +  E 


Once  again,  as  in  the  case  of  strain  hardening  behavior,  increments  of  displacement  and  total 
strain  are  obtained  from  the  linear  matrix  equations.  The  solution  for  from  Equation  27 
and  its  substitution  into  Equations  25  and  26  represent  the  procedure  necessary  to  obtain  the 
complete  solution  for  displacements,  stresses,  and  strains  for  any  increment,  assuming 
elastic,  ideally  plastic  material  behavior. 
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SECTION  m 

MATERIAL  AND  GEOMETRIC  NONLINEARITY 

For  the  preceding  theory  and  applications,  it  was  assximed  that  the  strain -displacement 
relations  were  linear  althou^  the  stress-strain  law  was  not.  The  implementation  of  a 
nonlinear  stress-strain  relation  for  the  characterization  of  material  behavior  merely  de¬ 
pends  on  the  absolute  magnitudes  of  the  elongations  and  shears  existing  in  a  body.  When 
they  exceed  a  certain  value,  nonlinear  material  characteristics  become  important  and  must 
be  Included  to  gain  an  insight  into  the  response  of  the  structure  to  further  loading,  Althou^ 
the  magnitudes  of  the  shears  and  elongations  may  be  sufficiently  large  to  necessitate  the 
inclusion  of  plastic  effects,  their  values  and  the  value  of  the  angles  of  rotation  may  still  be 
small  con^ared  to  unity.  If  this  condition,  and  the  additional  condition  that  the  scpiares  and 
products  of  angles  of  rotation  may  be  neglected  as  compared  to  the  elongations  and  shears 
remain  valid,  then  the  use  of  linear  strain -displacement  relations  is  justified.  Thus,  material 
nonlinearity  can  exist  independently  of  geometric  nonlinearity. 

For  flexible  bodies  (beams,  plates,  shells),  the  second  condition  (that  on  the  squares  and 
products  of  rotation)  is  not  satisfied  in  many  Instances.  Under  these  circumstances  it  is 
imjustiflable  to  neglect  the  terms  containing  the  squares  and  products  of  the  rotations  In  the 
strain-displacement  relations.  Furthermore,  the  linear  equilibrium  equations  are  no  longer 
valid  and  nonlinear  terms  consistent  with  the  Inclusionof  rotations  in  the  strain-displacement 
relations  should  be  retained.  Thus,  stresses  that  multiply  rotations  should  not  automatically 
be  dropped  In  deference  to  those  that  appear  linearly  in  the  equilibrium  equations. 

Geometric  nonlinearity  can  exist  independently  of  physical  nonlinearity  since  small  shears 
and  elongations  do  not  imply  small  angles  of  rotation.  Problems  requiring  the  consideration 
of  geometric  nonlinearity  alone  include  the  question  of  stability  of  elastic  equilibrium,  the 
deformation  of  bodies  having  initial  stresses,  large  deflection  of  beams,  plates  and  shells, 
and  torsion  and  bending  in  the  presence  of  axial  forces.  For  these  situations,  the  effect  of 
geometric  nonlinearity  must  be  taken  into  account  not  only  In  the  strain  displacement  relations, 
but  in  determining  changes  in  the  length  of  line  elements,  and  in  formulating  the  conditions 
of  equilibrium  of  the  volume  element.  In  addition,  if  the  magnitudes  of  the  strains  become  too 
large,  it  then  becomes  necessary  to  include  material  nonlinearity  through  the  stress-strain 
relations. 
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In  the  following,  the  procedure  developed  for  material  nonlinearity  is  expanded  to  include 
geometric  nonlinearity.  Althou^  only  elastic  perfectly-plastic  results  are  given,  the  method 
is  equally  applicable  and  easily  adaptable  to  strain-hardening  behavior  using  the  procedure 
outlined  In  the  previous  section.  Representative  structures  chosen  to  illustrate  the  signif¬ 
icant  features  of  combined  geometric  and  material  nonlinearity  are  the  restrained  beam  and 
the  circular  arch. 

Martin  (Reference  6)  has  presented  an  Incremental  numerical  method,  based  on  the 
direct  stif&iess  approach,  which  is  generally  applicable  for  the  treatment  of  problems  in¬ 
volving  geometric  nonlinearity.  This  procedure  approximates  the  nonlinear  behavior  by  a 
sequence  of  linear  steps.  Either  loading  or  displacement  may  be  applied  incrementally.  This 
procedure  also  requires  the  introduction  of  an  Initial  stress  stiffness  matrix,  in  addition 
to  the  conventional  stiffness  matrix,  to  account  for  the  effects  on  the  bending  stiffiiess  of 
the  membrane  loads  generated  as  a  consequence  of  geometric  nonlinearity;  that  is,  the  effects 
of  rotations  on  strains.  Thus,  the  implementation  of  this  matrix  in  addition  to  the  initial 
strain  matrix  represents  the  required  modifications  for  the  development  of  an  incremental 
procedure  to  account  for  both  types  of  nonlinearity. 

A  general  development  of  the  matrices  discussed  above  and  the  method  of  solution  to  the 
problem  of  combined  nonlinearity  is  obtained  by  following  the  procedure  as  outlined  in  Ref¬ 
erence  6  with  modifications  associated  with  the  inclusion  of  plasticity.  The  total  elastic  strain 
may  be  written  as  the  sum  of  three  components;  i.e.. 


•*  =  #°  +  A#^  -  A« 


(28) 


where  is  the  initial  elastic  strain  vector  (equal  to  P  /AE  for  the  beam  column,  and  re- 

o  ^ 


lated  to  the  initial  stresses  cr^,  a 


lated  to  the  initial  stresses  cr^,  cr^,  for  the  two  dimensional  plate  problem).  The  vector 
rp  X  y  xy 


rp  A  jr  Ajr 

Ae  is  the  additional  total  strain  developed  within  the  increment  of  load.  This  strain  in- 


T 

Ae  is  the  additional  total 


crement  is  related  to  the  Increments  of  displacement  throu^  the  strain-displacement  relations 
which  must  now  include  the  nonlinear  rotation  terms.  For  a  beam  column  element,  this  re¬ 
lationship  is  given  by 


(29) 
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where  Au  and  Aw  represent  the  increments  in  the  axial  and  lateral  displacements  of  the 
middle  surface  of  the  beam.  The  usual  Bernoulli -Euler  kinematic  beam  theory  assumptions 
were  made  to  obtain  Equation  29.  The  first  term  in  the  above  equation  represents  the  ex¬ 
tension  of  the  centerline  of  the  beam,  the  second  term  is  the  contribution  to  the  extensional 
strain  due  to  lateral  deflection  (the  rotation  term),  and  the  last  term  is  the  conventional 
bending  strain  term  arising  from  the  condition  that  normals  to  the  neutral  axis  should 
remain  after  deformation  straight  and  normal  to  the  centerline  and  unextended.  The  corre¬ 
sponding  strain-displacement  relations  for  a  plate  are 


Ae 


T 

X 


Ae 


T 

y 


_  d(Au) 
ax 

.  1  /'a(Aw)\^  a^(Aw) 

2  "  a«2 

a  ( Av) 

1  /dlAw)\  a^(Aw) 

dy 

+  2  [-^)  -  ^ 

aiAu) 

d(Av)  d  (Aw)  d(Aw) 

dy 

dx  dx  dy 

2 

d  (Aw) 

Sk  dy 


(30) 


The  functional  form  for  the  increments  in  displacement  for  the  beam  column  Is  chosen 
to  be  identical  to  that  given  In  Equations  1  and  4  for  the  total  displacement;  that  is, 

HqI*  (x)Aw.  +  h|V  (x)Aw,^ 

i  _ 

2  (3i; 

Au(x)  =  z  Hq?^(x)Au 
i  =  l 

where  the  definitions  of  the  Hermitian  polynomials  may  be  obtained  by  comparing  corre¬ 
sponding  terms  in  Equations  1,  4,  and  31.  A  representation  identical  to  that  given  in  Ref¬ 
erence  9  may  be  used  for  the  increment  of  lateral  displacement  for  the  plate  element.  The 
incremental  relationship  between  total  strain  and  generalized  nodal  displacements  is  obtained 
by  taking  appropriate  derivatives  of  the  assumed  displacement  functions. 


Aw(x)  =  Z 

i  =  I 


The  increment  in  elastic  strain  energy  from  an  initial  elastic  strain  state,  may  be 
written  as 


AU  =  /// 


/ 


•  °+  A«® 


<r' 


dV 


(32) 
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Integrating  the  above  equation  between  the  prescribed  limits  of  strain  and  using  the  following 
relations: 


o-  =  E 
cr®  =  E 


and 

A«*  =  A#  ^  -  A<, 

we  can  write  the  increment  of  strain  energy  as 

T  *  T  T* 

E  A#  dV  -  ///A.'  E  Ac  dV  + 

''  (331 

;//.«■  A/  dV  +  i-  Iff  A.'  E  A.  dV  -  ///.»■  A.  dV 

V  V  V 

T 

Only  those  terms  In  the  above  equation  that  contain  ,  which  is  a  function  of  the  Increment 
of  nodal  displacement,  will  contribute  to  the  incremental  load-deflection-inltial  strain 
relationship.  The  remaining  terms  are  arbitrary  constants.  Upon  neglecting  higher  order 
terms  in  the  increments  of  displacement,  the  first  term  in  Equation  33  leads  to  the  con¬ 
ventional  stiffness  matrix.  The  second  term  yields  the  initial  strain  matrix,  while  the  third 
term  leads  to  the  initial  stress  stiffness  matrix.  In  the  development  of  the  latter  matrix  the 
work  done  by  the  in -plane  stresses  and  the  generation  of  additional  membrane  stresses, 
resulting  from  the  effects  of  geometric  nonlinearity,  are  both  taken  to  be  zero  during  the 
application  of  an  Increment  of  lateral  load.  These  considerations  constitute  the  linearization 
of  the  procedure  during  an  increment  of  load. 

Consistent  with  the  incremental  procedure  of  the  previous  section,  the  increment  of 
plastic  strain  may  be  written  as 


Ac’ 


C  -  € 


I  -» 


Substituting  the  above  equation  into  Equation  33  and  then  making  use  of  Casti^iano’s  theorem 
leads  to  the  following  relation  for  an  individual  element: 


ARn 


JOJ  A.  > 

k  Ado 


+  k''»Ado' 


-  k 


i-  I 


(34) 
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where 

=  fff  W‘  E  W  dV, 
V 


the  conventional  stiffness  matrix, 

V 

the  initial  stress  stiffness  matrix,  and 

k*  '  ///W  E  Wp  dV  . 

the  initial  strain  stiffness  matrix.  The  matrix  W  relates  the  increment  of  rotation  to  the 
Increment  of  generalized  nodal  displacement  and  therefore  represents  the  nonlinear  con¬ 
tribution  to  the  strain-displacement  relations.  The  matrices  W  and  Wp  are  the  same  as  those 
used  in  the  previous  section, 

A  predictor-procedure  must  once  again  be  used  to  obtain  a  solution  because  the  location 
of  the  elastic-plastic  boundary  is  not  known  a  priori.  Thus,  the  governing  linear  matrix 
equation  is  written  in  the  following  form: 

/,(0)  ,.(1)  \  .  i  .  i  .1-1 

U  ^  I  <25) 

where  Aq  is  defined  in  Equatipn  13,  and  must  once  again  be  retarded  by  one  step  in  the 
solution  procedure.  Equation  35  is  very  nearly  the  same  as  Equation  13  used  in  Section  II 
for  the  plastic  bending  analysis.  However,  because  of  the  presence  of  geometric  nonlinearity, 
the  entire  element  stiffness  matrix  k  =  must  be  reformed  at  every  step  using 

current  stress  levels  and  geometry. 

Thus  the  solution  procedure  requires  that  for  a  generic  step,  and  are  calculated 
by  making  use  of  the  geometry  and  initial  forces  existing  at  the  start  of  the  step.  The  in¬ 
crement  in  the  fictitious  force,  Aq  ,  is  calculated  using  current  and  the  immediately  preceding 
values  of  the  location  of  the  elastic-plastic  boundary.  An  Increment  of  load  is  then  applied, 
and  the  corresponding  displacement  increments  calculated  from  the  total  matrix  equation 
obtained  by  assembling  Equation  35  for  each  element.  New  internal  forces  are  calculated 
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and  total  stresses,  strains,  and  displacements  are  obtained  by  summing  Incremental  values. 
The  new  location  of  the  elastlc-'plastlc  boundary  is  determined  from  the  total  strain  (or  stress) 
distribution  and  the  process  is  repeated  until  the  maximum  specified  load  level  is  reached  or 
the  structure  falls. 


SECTION  IV 
RESULTS 

As  a  demonstration  of  the  feasibility  of  the  plastic  bending  analysis,  application  of  the 
method  has  been  made  to  some  elementary,  but  representative  sample  structures.  For  two  of 
these  structures  (a  simply  si^jported  and  a  cantlleveredbeam),  results  from  an  exact  solution 
to  the  governing  differential  equation,  assuming  elastic-ldeally  plastic  material  behavior,  are 
available  for  comparison.  As  a  consequence  of  assuming  elastic-ldeally  plastic  behavior, 
and  since  both  of  the  structures  are  statically  determinate,  an  analytic  expression  can  be 
written  which  relates  the  depth  of  the  elastic -plastic  boundary  to  the  sqjplied  load.  The  finite 
element  analysis  is  initially  applied  to  the  beam  structures  using  this  relationship,  thus  pro¬ 
viding  a  means  of  determining  the  validity  of  assumptions  made  in  choosing  such  quantities  as 
the  displacement  function,  the  plastic  strain  distribution,  and  the  representation  of  the 
elastic-plastic  boundary. 

Figure  4a  represents  a  nondimensionallzed  load  versus  central  deflection  curve  for  a 
uniformly  loaded,  simply  supported  beam.  Six  elements  are  used  in  the  idealization  of  one- 
half  of  the  beam.  In  this  figure,  w^  is  the  center  deflection,  w^*  is  the  center  deflection  at 
the  maximum  load  for  which  the  beam  is  entirely  elastic,  and  represents  the  nondimenslonal 
load  parameter,  given  as 
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Here,  p  is  the  applied  load,  and  p^  =  ^^^yield’  results  obtained  from  the  finite  element 
analysis  compare  quite  favorably  with  corresponding  results  from  the  exact  solution  (Ref¬ 
erence  10),  as  shown  in  Figure  4a,  The  collapse  load,  as  determined  from  the  near  vertical 
slope  of  the  load-deflection  curve,  is  approximately  three  percent  higher  than  the  exact 
collapse  load  which  occurs  at  a  value  of  p  -  1. 


The  progression  of  the  elastic -plastic  boundary  through  the  thickness  and  in  the  plane  of 
the  element  is  shown  in  Figure  4b,  From  this  figure  It  is  seen  that,  although  the  depth  of  the 
boundary  at  plastic  nodes  is  exact,  the  assumption  associated  with  its  distribution  in  the  plane 
(i,  e.,  linearly  varying  to  adjacent  nodes),  may  lead  to  discontinuities  In  the  representation  of  its 
slope  as  evidenced  in  the  figure  for  load  values  oi  p  =  1.00  and  P-  1,03.  The  appearance  of 
these  discontinuities  indicates  that  the  actual  boundary  lies  between  the  nodes  on  the  upper  and 
lower  surface  of  the  beam.  The  error  Introduced  by  the  assumption  of  a  linearly  varying 
boundary  in  the  plane  can  be  reduced  by  increasing  the  number  of  elements  used  In  the 
idealization  of  the  beam.  Also  to  be  noted  in  Figure  4b  is  the  development  of  a  fully  plastic 
cross  section  located  at  the  center  of  the  beam  at  a  load  of  =  1,  In  a  continuum  analysis, 
the  development  of  this  fully  plastic  cross  section  is  sufficient  to  cause  collapse  of  this 
structure.  However,  in  the  finite  element  analysis,  coll£q)Be  does  not  occur  until  both  cross 
sections  of  the  element  containing  the  center  of  the  beam  become  fully  plastic. 

Results  in  the  form  of  a  nondimensionallzed  load  versus  tip-deflection  curve  for  a  uni¬ 
formly  loaded,  cantilevered  beam,  are  shown  In  Figure  5a.  Elastic-ideally  plastic  material 
behavior  was  assumed.  A  comparison  with  results  from  an  exact  solution,  shown  as  the  solid 
curve  in  the  figure,  indicates  good  correlation  up  to  the  collapse  load.  For  this  problem,  as 
for  the  simply  supported  beam,  the  depth  of  the  elastic-plastic  boundary  can  be  directly 
related  to  the  ^plied  load.  This  relationship  was  used,  once  again,  to  obtain  the  results  shown 
in  the  figure. 

The  progression  of  the  elastic-plastic  boundary,  throu^  the  cross  section  and  in  the  plane 
of  the  elements,  is  shown  in  Figure  5b.  As  indicated  in  the  figure,  the  development  of  the 
boundary  is  much  more  localized  for  this  structure  than  it  was  for  the  sln^)ly  supported  beam. 
Consistent  with  a  continuum  approach,  collapse  of  this  structure  was  indicated  in  the  finite 
element  analysis  by  the  development  of  one  fully  plastic  cross  section  at  the  root  of  the 
cantilever. 
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Figure  4a.  Load  Versus  Deflection  of  a  Simply  Supported  Beam 


Figure  4b.  Progression  of  Elastic-Plastic  Boundary:  Simply  Supported  Beam 
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Figure  5a.  Load  Versus  Deflection  of  a  Cantilevered  Beam 


Figure  5b.  Progression  of  Elastic-Plastic  Boundary  Cantilever  Beam 
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For  both  the  simply  supported  beam  and  the  cantilevered  beam,  as  previously  mentioned, 
an  exact  relationship  between  the  depth  of  the  elastic -plastio  boundary  at  nodes  in  the  plastic 
range  and  the  applied  load  was  used  to  obtain  results  using  the  finite  element  analysis.  The 
Justification  for  using  this  relationship,  which  admittedly  does  not  exist  for  most  structures 
of  interest,  was  to  check  the  validity  of  assumptions  made  in  choosing  such  quantities  as  the 
displacement  function,  the  plastic  strain  distribution,  and  the  representation  of  the  elastic- 
plastic  boundary.  As  Indicated  by  the  previous  results,  the  use  of  these  assumptions  for  the 
finite  element  analysis  appears  to  be  justified. 

Since,  In  general,  the  depth  of  the  elastic-plastic  boundary  Is  not  known  at  the  current 
load  step,  results  for  the  cantilevered  beam  were  recomputed  and  a  load-deflectlon  curve 
obtained,  by  using  an  approximate  value  for  the  depth  of  the  elastic -plastic  boundary  which  Is 
shown  in  Figure  6.  The  value  of  the  depth  of  the  boundary  used  for  any  increment  of  load  Is 
based  on  the  total  strain  distribution  determined  at  the  end  of  the  previous  load  Increment. 
The  use  of  this  procedure  cannot  lead  to  the  development  of  a  fully  plastic  cross  section. 
Consequently,  It  Is  assumed  that  a  fully  plastic  cross  section  exists  at  a  node  when  plasticity 
has  developed  a  specified  amount  through  the  thickness.  The  deflections  and  the  slope  of  the 
load -deflection  curve  for  this  structure,  increase  quite  rapidly  beyond  a  value  of  load  for 
which  plasticity  has  developed  through  80  percent  of  the  end  cross  section.  Thus,  In  the  analysis 
this  value  was  chosen  as  the  criterion  to  determine  the  development  of  a  fully  plastic  cross 
section.  The  degree  of  approximation  obtained  by  using  previous  values  of  the  depth  of  the 
boundary,  when  compared  with  the  exact  solution.  Is  Illustrated  in  Figure  6.  As  can  be  seen, 
the  results  compare  favorably  for  most  of  the  load  range  considered.  The  maximum  diver¬ 
gence  occurs  in  the  vicinity  of  the  collapse  load,  determined  by  using  the  exact  boundary 
position.  The  error  introduced  by  using  the  approximate  position,  for  the  prediction  of  the 
collapse  load.  Is  seven  percent. 

Results  are  also  shown  in  Figure  6  for  the  cantilevered  beam,  for  which  strain  hardening 
material  behavior  Is  assumed.  These  results,  shown  as  the  dotted  curve,  are  compared  with 
the  corresponding  results  obtained  by  using  elastic -Ideally  plastic  behavior.  The  proximity 
of  results  for  strain  hardening  and  perfectly  plastic  behavior  can  be  attributed  to  the  use  of 
Ramberg-Osgood  strain  hardening  parameters,  chosen  to  approximate  the  elastic-ldeally 
plastic  stress-strain  curve.  The  slope  of  the  load-deflection  curve  for  strain  hardening 
behavior  illustrates  that  the  beam  still  possesses  some  stiffness  beyond  the  theoretical 
collapse  load  predicted  by  assuming  perfectly  plastic  behavior. 
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Figure  7  illustrates  the  application  of  the  procedure  to  a  simply  supported  beam  sub¬ 
jected  to  combined  bending  and  axial  loads.  As  previously  discussed  in  Section  III,  the  analysis 
for  this  problem  requires  the  introduction  of  an  Initial  stress  stiffness  matrix  to  account  for 
the  effects  of  the  membrane  load  on  the  bending  stiffness.  The  determination  of  the  position  of 
the  elastic -plastic  boundary  relative  to  both  the  upper  and  lower  surfaces  Is  also  required 
for  this  analysis.  Results  have  been  obtained  for  cases  in  which  a  uniform  lateral  load  acts 
In  conjunction  with  a  constant  tensile  or  compressive  axial  load,  Indicated  in  the  figure  by 
T  =  +1000  and  T  ='  -1000,  respectively.  These  results  are  compared  with  the  case  of  pure 
bending.  Indicated  as  T  =*  0,  As  seen  in  the  Figure  7,  the  effect  of  the  axial  compressive  load 
is  to  reduce  the  stiffness  of  the  structure,  and  the  tensile  load  increases  the  stiffness,  when 
compared  to  the  case  of  pure  bending.  A  solution  to  this  problem  by  using  a  continuum 
analysis  similar  to  the  one  developed  for  pure  bending  In  the  plastic  range  does  not  appear 
to  be  available  for  comparison.  For  the  case  of  the  compressive  axial  load,  the  lateral  load 
was  Incremented  to  a  value  that  resulted  in  the  failure  of  the  structure.  This  failure  is  indi¬ 
cated  in  Figure  7  by  the  near  vertical  slope  of  the  load-deflection  curve.  It  should  be  noted 
that,  for  this  problem,  it  was  not  necessary  to  develop  a  completely  plastic  cross  section  for 
collapse  to  occur.  The  reduction  of  the  stiffness,  caused  by  the  axial  compressive  load  and 
the  progression  of  the  elastic-plastic  boundary  through  only  a  portion  of  the  thickness,  was 
sufficient  to  cause  failure.  This  type  of  failure  is  associated  with  plastic  buckling  rather  than 
the  formation  of  a  mechanism. 

Application  of  this  procedure  has  been  made  to  a  simply  supported,  uniformly  loaded, 
square  plate.  Using  a  36-element  Idealization  to  represent  the  quarter  panel,  load  versus 
central  deflection  curves  for  this  structure,  assuming  elastic-ldeally  plastic  and  elastic- 
strain  hardening  material  behavior,  are  shown  In  Figure  8.  Once  again,  as  in  the  case  of  the 
beam,  the  proximity  of  results  for  both  types  of  material  behavior  is  attributable  to  the  choice 
of  strain  hardening  parameters  that  approximate  elastic-ldeally  plastic  material  behavior. 

The  collapse  load  for  this  structure,  determined  by  assuming  elastic-ldeally  plastic 
material  behavior.  Is  the  value  of  the  load  at  which  the  pattern  of  fully  plastic  elements  forms 
a  mechanism  of  failure.  The  pattern  of  the  development  of  the  plastic  region  in  the  plane  of 
the  plate,  and  the  progression  of  the  boundary  through  the  thickness  of  the  plate.  Is  shown  In 
Figures  9a  and  9b,  respectively.  In  Figure  9a  the  cross-hatched  area  represents  those 
regions  of  the  plate  In  which  plasticity  has  developed  to  some  degree  less  than  80  percent 
through  the  cross  section.  The  shaded  area  represents  those  regions  in  which  plasticity  has 
developed  to  a  depth  greater  than  80  percent  throu^  the  thickness.  A  consideration  of  the 
latter  region  as  being  fully  plastic  leads  to  the  development  of  a  mechanism  of  collapse 
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Figure  6.  Load  Versus  Tip  Deflection  of  Cantilever  Beam 


Figure  7.  Simply  Supported  Beam  Subjected  to  Combined  Bending  and  Axial  Loads 
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Figure  8.  Load  Versus  Central  Deflection  of  a  Uniformly  Loaded  Simply  Supported 
Square  Plate 
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formed  along  the  diagonals  of  the  square  plate,  as  shown  In  Figure  9a,  As  In  the  case  of  the 
beam,  this  criterion  Is  necessary  because  the  depth  of  the  boundary,  as  cort^Juted  from  the 
total  strains,  cannot  lead  to  the  development  of  a  fully  plastic  section. 

The  pattern  of  development  of  the  plastic  region  In  the  plane  of  a  narrow  rectangular 
plate  -  0.3)  and  the  progression  of  the  boundary  through  the  thickness  across  the  length 
and  width  of  the  plate,  is  shown  in  Figures  lOa  throu^  10c,  respectively.  In  Figure  10a  the 
80  percent  criterion  was  once  again  used  to  determine  the  pattern  of  fully  plastic  sections 
in  forming  the  collapse  mechanism.  From  this  figure  It  Is  evident  that  the  sections  along  the 
collapse  pattern  do  not  all  lie  on  the  diagonals  of  the  plate. 

A  comparison  of  upper  bound  solutions  for  the  load  carrying  capacities  of  rectangular 
plates  of  various  aspect  ratios  is  shown  in  Figure  11.  The  solid  curve  represents  the  solution 
obtained  using  the  von  Mises  yield  criterion  in  conjunction  with  assumed  collapse  pattern  (1), 
shown  In  Figure  11.  The  dotted  curve,  obtained  from  Reference  12,  represents  the  t^per 
bound  solution  obtained  using  the  Tresca  yield  condition,  together  with  assumed  collapse 
pattern  (2),  Results  from  the  finite  element  analysis  are  represented  by  the  solid  circles. 
The  finite  element  results  indicate  that  the  displacement  pattern  (2)  Is  a  more  accurate 
representation  of  the  collapse  mechanism  than  is  pattern  (1),  An  upper  bound  solution  using 
the  second  displacement  pattern  with  the  von  Mises  yield  condition  is  shown  as  the  dashed 
curve  in  Figure  11.  The  results  from  the  present  analysis  compare  favorably  with  this  latter 
solution  and  are  slightly  below  It  except  for  extremely  low  aspect  ratios.  For  such  narrow 
plates,  the  use  of  the  80  percent  criterion  in  conjunction  with  the  depth  of  the  elastic -plastic 
boundary,  calculated  from  the  total  strain  distribution  of  the  previous  step,  is  not  adequate, 
A  relaxation  of  the  80  percent  criterion,  based  on  a  careful  examination  of  the  load  deflection 
history,  spears  to  be  warranted.  In  addition,  a  possible  alternative  might  be  the  incorporation 
of  an  Iterative  procedure  in  the  method  of  solution. 

To  illustrate  the  procedure  associated  with  geometric  nonlinearity  and  combined  material 
and  geometric  nonlinearity,  a  simply  supported,  restrained  beam,  subjected  to  a  uniform 
vertical  load,  is  considered.  Load  versus  central  deflection  curves,  obtained  for  purely 
elastic  and  for  elastic -ideally  plastic  material  behavior  are  shown  in  Figure  12a,  Results 
assuming  elastic  behavior  from  the  finite  element  analysis  and  exact  results  from  Ref¬ 
erence  13  are  shown  as  the  dotted  curve  in  the  figure.  Load  versus  deflection  curves,  as¬ 
suming  plastic  behavior  are  presented  for  idealizations  involving  6,  12,  and  24  elements  for 
one-half  of  the  beam.  Differences  in  the  results  for  these  idealizations  only  appear  after  the 
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Figure  9.  Progression  of  Elastic-Plastic  Region  for  Simply-Supported,  Uniformly  Loaded, 
Square  Plate 
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Figure  10.  Progression  of  Elastic-Plastic  Region  for  Simply-Supported,  Uniformly 
Loaded  Rectangular  Plate 
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end  sections  at  the  supports  become  fully  plastic.  Deflections  beyond  the  value  of  load  at 
which  this  occurs  increase  quite  rapidly,  and  collapse  occurs  shortly  thereafter  with  the 
development  of  another  fully  plastic  cross  section.  The  counterbalancing  effect  of  geometric 
and  material  nonlinearity  is  vividly  depicted  In  Figure  12a  by  the  region  of  the  load-deflection 
curve  which  is  very  nearly  linear. 

Figure  12b  illustrates  the  growth  of  the  plastic  regions  of  the  restrained  beam.  The 
dotted  line  at  P  ==  179  Ib/ln.  indicates  the  jump  in  the  representation  of  the  plastic  region, 
when  the  end  section  becomes  fully  plastic. 

The  load-deflection  history  of  a  circular  arch  subjected  to  a  concentrated  load  is  shown 
in  Figure  13.  The  elastic  buckling  load  compares  favorably  with  that  obtained  in  Reference  6. 
Load  versus  center  deflection  curves,  obtained  by  assiuning  elastic-ideally  plastic  material 
behavior,  are  shown  for  two  values  of  yield  stress.  The  onset  of  collapse  for  this  structure  is 
appreciably  hastened  with  the  introduction  of  plasticity.  This  is  attributable  to  the  reduction 
of  the  stiffiiess  of  the  structure  resulting  from  the  complementary  effects  of  geometric  and 
physical  nonlinearity.  As  in  the  case  of  the  uniformly  loaded  beam  subjected  to  a  constant 
axial  compressive  load,  the  development  of  a  fully  plastic  cross  section  was  not  necessary 
for  collapse  to  occur. 

Figure  14  Illustrates  the  load-deflection  histories  of  the  same  arch  as  that  used  in 
Figure  13,  now  subjected  to  a  uniform  load  distribution.  Once  again  it  is  seen  that  the  effect 
of  plasticity  considerably  reduces  the  collapse  load  of  the  structure  from  its  elastic  buckling 
load. 
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Figure  12a,  Elastic-Plastic  Response  of  a  Restrained  Beam 


I 


Figure  12b,  Progression  of  Elastic-Plastic  Boundary  In  Restrained  Beam  for  Increasing 
Load 
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Figure  13.  Load  Versus  Central  Deflection  for  a  Simply  Supported  Arch  Subjected  to  a 
Concentrated  Load 


Figure  14.  Load  Versus  Central  Deflection  of  a  Sinyply  Supported  Arch  Subjected  to  a 
Uniform  Load 
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SECTION  V 

SUMMARY  AND  CONCLUSIONS 

A  finite  element  method  Is  presented  that  can  account  for  material  nonlinearity,  alone, 
or  in  combination  with  geometrically  nonlinear  behavior  for  the  out-of-plane  bending  analysis 
of  structures.  The  initial  strain  concept  is  implemented  into  the  finite  element  analysis, 
formulated  within  the  framework  of  the  direct  stiffaess  method,  to  accoxmt  for  the  effects  of 
plasticity.  These  effects  are  introduced  into  the  analysis  as  a  fictitious  load  vector  to  be 
combined  with  the  actual  applied  load.  Thus,  the  method  of  analysis  can  be  readily  incor¬ 
porated  into  existing  finite  element  procedures. 

Application  of  the  method,  illustrating  the  plastic  behavior  of  typical  structures  under 
pure  bending,  and  bending  in  combination  with  applied  axial  loads,  is  presented.  Results 
from  the  finite  element  analysis  are  compared  with  results  from  analytical  solutions, 
where  possible,  for  beam  and  plate  structures.  Good  correlation  is  indicated  for  the  load- 
deflection  characteristics  of  these  structures  as  well  as  for  the  prediction  of  plastic  collapse 
loads.  In  addition,  with  the  generation,  or  application  of  compressive  membrane  stresses, 
the  present  procedure  is  capable  of  predicting  failure  resulting  from  a  combination  of  plastic 
collapse  and  buckling.  The  correlation  of  results  and  the  numerical  stability  of  the  procedure 
as  applied  to  the  sample  problems  substantiates  the  assumptions  associated  with  the  func¬ 
tional  representation  of  the  elastic -plastic  boundary  and  the  distribution  of  the  plastic  strain 
within  each  element,  in  addition  to  the  use  of  the  predictor  form  of  the  solution  procedure. 

Since  the  phenomenon  of  plastic  deformations  may  itself  lead  to  large  displacements 
and  rotations,  the  treatment  of  effects  arising  from  geometric  nonlinearity  assumes  partic¬ 
ular  significance  in  the  solution  of  many  important  problems.  Consequently,  the  plastic 
bending  analysis  was  combined  with  a  method  capable  of  accounting  for  geometrically  non¬ 
linear  behavior.  This  combined  procedure  was  applied  to  a  restrained  beam  and  a  simply 
supported  circular  arch.  In  the  beam,  the  effects  of  geometric  nonlinearity  act  counter  to 
the  reduction  of  stiffness  caused  by  the  progressive  development  of  plasticity.  Failure  of  this 
structure  occurs  only  after  the  development  of  fully  plastic  cross  sections.  For  the  circular 
arch,  the  effects  of  both  types  of  nonlinearity  are  complementary,  and  failure  occurs  as  a 
plastic  buckling  phenomenon. 
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APPENDIX 

INITIAL  STRAIN  STIFFNESS  MATRICES 
FOR  BEAM  FINITE  ELEMENTS 

The  initial  strain  stiffness  matrices  for  abeam  element  in  pure  bending  and  for  combined 
bending  and  membrane  loading,  are  derived  on  the  basis  of  the  assumptions  shown  in  Figures  1 
and  2,  and  are  given  in  integral  form  in  Equation  10. 


The  matrix  equation  defining  the  fictitious  nodal  restoring  forces  in  terms  of  the  re¬ 
sulting  initial  strain  stiffness  matrix,  is  shown  below  for  the  pure  bending  of  a  beam  with  a 
rectangular  cross  section,  represented  as 
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- ^ - - —  + - i - - - ! - +—! - !- 


20 


12 
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and  P  and  M  represent  the  fictitious  restoring  force  in  the  lateral  direction  and  moments, 
z 

respectively.  For  this  element,  all  T’s  are  determined  with  respect  to  the  median  surface. 
Other  quantities  appearing  in  Equation  36  are  defined  in  Figure  1. 


The  corresponding  relations  for  the  case  of  combined  bending  and  membrane  stresses 
are  shown  below: 


(37) 


where 

k*  =  -  {jt^/40)c,  +  (X/4)c3 

k*2  =  -  (9jt^/40)c,  -  (1^/4)  Cg  -  (i/4)  Cj 

k*  =  -  {i*/40)c  +  (//4)c^ 

k,*  =  -  (9i’/40)c^  -  (i^/4)  Cj  -  (i/4)Cg 

k*,  =  {ji^/l20)c,  +  {i^/24)C2  +{/^/4)Cj 

kjg  =  -  {i^/20)c,  -  (/^/24)C2 

kgj  =  {Jt'*/I20)c^  +  (Jt’/24)Cg -f  (i^/4)Cg 

k*^  =  -  {i^20)c^  -  (i^/24)Cg 

k  *  =  -  (Z^/  8)  -  {3i^/8)Ca 
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k ‘j  =  -  (t’z  A)  c  j  -  (3/^/  8>Cj 

S*  =  +  I3^'/8IC,„ 

+  (3/"/8lc„ 

k*  =  -  (£‘’/30)c  -  (/’/24)c 

51  I  Z 

k*2  =  -  (7£‘’/40)c,  -  (5i^/24)C2  -(A^/4)Cj 

k*j  =  -  (je^/30)c^  -  (jt^/24)Cg 

k*^  =  -  (7XV40)C^  -  (5X*/24)Cg  -  (i^lCg 


and 
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c 


8 


i 


c 


9 


// 


c 


10 


-  2t 


The  terms,  P  and  P  ,  represent  the  fictitioiis  restoring  forces  in  the  axial  direction. 

Xi 

Here  all  z’s  are  measured  with  respect  to  the  upper  surface  of  the  beam.  All  other  quantities 
appearing  in  Equation  37  are  defined  In  Figure  2, 
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